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g : Abstract 
C — . 

f-*) \ We present an ansatz for the planar five-loop four-point amplitude in maximally supersym- 



metric Yang-Mills theory in terms of loop integrals. This ansatz exploits the recently observed 
correspondence between integrals with simple conformal properties and those found in the four- 
point amplitudes of the theory through four loops. We explain how to identify all such integrals 
systematically. We make use of generalized unitarity in both four and D dimensions to determine 
the coefficients of each of these integrals in the amplitude. Maximal cuts, in which we cut all 
propagators of a given integral, are an especially effective means for determining these coefficients. 
The set of integrals and coefficients determined here will be useful for computing the five-loop cusp 
anomalous dimension of the theory which is of interest for non-trivial checks of the AdS / CFT dual- 
ity conjecture. It will also be useful for checking a conjecture that the amplitudes have an iterative 
structure allowing for their all-loop resummation, whose link to a recent string-side computation 
by Alday and Maldacena opens a new venue for quantitative AdS/CFT comparisons. 
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I. INTRODUCTION 



Maximally supersymmetric Yang- Mills theory (MSYM) is an important arena for explor- 
■ng the properties of gauge theories. The Maldacena weak-strong duahty |1] between MSYM 

and string theory in AdSs x provides an explicit realization of 't Hooft's old dream [2(] 
of expressing the strongly coupled limit of a gauge theory in terms of a string theory. In 
addition, the higher-loop planar space-time scattering amplitudes of MSYM appear to have 
a remarkably simple and novel iterative structure 3, J]. This structure allows higher-loop 



amplitudes to be expressed in terms of lower-loop amplitudes. This simpjicity may well be 

LZl. 



connected to the observed integrability of the theory in the planar limit 

The iterative structure of the planar amplitudes was first proposed in ref. [3|], and con- 
firmed for the two-loop four-point amplitude. An independent two-loop check was given in 

n n 

ref. [8[. In ref. [4], this proposal was fleshed out for all planar maximally helicity violating 
(MHV) amplitudes, expressed via a specific all-loop exponentiation formula which was con- 
firmed for three-loop four-point amplitudes. The iteration formula has also been shown to 
hold for two-loop five-point amplitudes [9] . At four loops the amplitudes are known in terms 
of a set of integrals, but the integrals themselves have not yet been evaluated fully 10, lll |- 
In this paper, we provide the corresponding integral representation of the five-loop four- 
point planar amplitude, for use in future studies of its properties. In a very recent paper, 
Alday and Maldacena 12j have shown how to perform a string-side computation of the same 
gluon amplitudes in the strong-coupling limit. This opens new and exciting possibilities of 
quantitative checks of the AdS/CFT correspondence, going beyond anomalous dimensions 
to detailed dependence on kinematics. 

In addition to exhibiting an iterative structure, the scattering amplitudes provide new 
and nontrivial information on the AdS/CFT correspondence. Using considerations of inte- 
grability, an integral equation for the cusp (soft) anomalous dimension — valid to all loop 
orders 

first three loop orders 

it would hold beyond this. We now know that this original proposal requires modification 
because of the recent calculation of the four-loop cusp anomalous dimension from the in- 



was written ^own by Eden and Staudacher This equation agreed with the 
but its reliance on various assumptions cast doubt on whether 



frared singular terms of a four-loop MSYM amplitude 10, lUl • A remarkable new integral 



equation proposed by Beisert, Eden and Staudacher (BES) ^ is in agreement with this 
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calculation. Surprisingly, the first four loop orders of the planar cusp anomalous dimension 
contain sufficient information to test the AdS / CFT correspondence to the level of a few per- 
cent [10|], u sing the interpolating function technique introduced by Kotikov, Lipatov, and 



Velizhanin 



14| as well as Fade approximants. This provides an independent guess of the 



entire perturbative series lOl], matching the one generated by the BES integral equation. 
Detailed studies of the BES equation confirm that it has the proper behavior [l^ at 
strong coupling, giving a high degree of confidence in it. Nonetheless, further checks on the 
perturbative side would be quite valuable. The first such check requires a computation of 
the five-loop anomalous dimension, which — following the approach taken at four loops — 



requires an expression for the five-loop 



The unitarity method [17|, 
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20 



bur-point amplitude. 
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22l | provides a powerful method for computing 



gauge and gravity loop 



amplitudes and has played a central role in obtainin g M SYM loop 



amplitudes 3|, Ij, |9|, llO. llTl . llSl . |23| . l2J, |25|]. An important recent improvement 26| is the use 
of complex momenta 27| within the framework of generalized unitarity 22|. In particular, 
this allows one to define a non-zero massless three-point amplitude, which vanishes for real 
momenta. At one loop this enables an easy algebraic determination of the coefficient of any 
box integral appearing in the theory, because the cut conditions freeze the loop integrals 26|. 
Some of these ideas have also been apphed at two loops [2^;]. In this paper we apply these 
ideas to develop a maximal- cutting method for efficiently determining coefficients of higher- 
loop integrals. 

In the MSYM theory, a special set of cuts — the iterated two-particle cuts — give rise to 
the "rung rule" for systematically obtaining higher-loop integral representations of planar 
four-point amplitudes 23|, l24] • At two and three loops this rule generates all contributions 
appearing in the amplitudes. However, starting at four loops, new integrals arise which are 
not generated by the rung were computed explicitly using generalized 

unitarity by relying on a set of mild assumptions. 

These integrals, along with two others that do not appear in the amplitude, are pre- 
dicted by a procedure relying on a beautiful observation due to Drummond, Henn, Smirnov 
and Sokatchev (DHSS) [291]. These authors noticed that, through three loops, the massless 
integrals appearing in the planar four-point amplitude are in direct correspondence with 
conformally invariant integrals. This correspondence comes from replacing dimensional reg- 
ularization with an off-shell infrared regularization in four dimensions. We shall call the 
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dimensional regularized version of the conformal integrals pseudo-conformal, since dimen- 
sional regularization breaks conformal invariance. DHSS also gave simple rules for generating 
all such integrals via "dual diagrams". The direct evaluation of generalized unitarity cuts 



confirmed 



10|] at four loops, that only such pseudo-conformal integrals appear in the planar 



amplitude. In the present work, we will assume that this is also true at higher loops, be- 
ginning with five loops. This provides a basis set of integrals for the planar (leading-color) 
contributions to the five-loop amplitude. We then use the unitarity method to determine 
the coefficients of these integrals in the planar five-loop four-point amplitude as well as to 
provide consistency checks on the absence of other integrals. 

We make use of an additional observation: the cutting equations hold at the level of the 
integrands, prior to carrying out any loop integrals. Indeed they hold independently for each 
of the multiple solutions of the cutting equations. These properties are especially powerful 
when combined with the basis of pseudo-conformal integrals. The problem is then reduced 
to an algebraic problem of determining the coefficient of each integral. Remarkably, it turns 
out that after dividing by the tree amplitude, the coefficients are pure numbers taking on 
the values —1, 0, or 1. This property is already known to hold for the four-point amplitude 
through four loops [lO^, and here we confirm it through five loops. 

At one loop, complete dimensionally regularized amplitudes in the MSYM theory can be 
constructed u sing only four- dimensional helicity amplitudes, greatly simplifying their con- 



struction 
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181]. Unfortunately, no such theorem exists at higher loops. Any rigorous 



construction of amplitudes requires that D-dimensional momenta be used in the cuts. It is 
worth noting that if our assumption of a pseudo-conformal basis of integrals is correct, with 
dimension independent coefficients, then unitarity in four dimensions does suffice to deter- 
mine these amplitudes in all dimensions. Our partial checks of dimensional cuts provide 
non-trivial evidence that this assumption is correct. This is rather remarkable because away 
from four dimensions there is a priori no reason why a simple analytic continuation of the 
dimension of the integrals should give correct results. 

Our expression for the five-loop four-point planar MSYM amplitude in terms of integral 
functions should be useful in a number of studies. The infrared singularities present in 
the amplitude encode the so-called cusp or soft anomalous dimension. At five loops these 
singularities begin at (where e, as usual, is the dimensional regularization parameter: 

e = (4 — D)/2). Evaluation of the infrared singular terms through would allow the 
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extraction of the five-loop cusp anomalous dimension, as has been done at three and four 
loops 0, ll^ . An evaluation though order 1/e would allow the extraction of a second anoma- 
lous dimension connected to a form factor. An evaluation through 0{e^) would allow a 
five-loop check of the iterative structure of the amplitudes, providing strong evidence that it 
continues to all loop orders. Although evaluating loop integrals is rather challenging, there 
has been rapid progress using Mellin-Barnes representations [sqI and in automating the re- 
quired analytic continuations jsil, allowing for explicit computations through four loops. 
Recently, there has also been progress in isolating the subsets of terms which determine 



the anomalous dimension 



greatly simplifying its calculation. Further development will 



34|. 



presumably be needed to apply these advances to the five-loop amplitude. 

Another important reason for studying MSYM amplitudes is their intimate connection 
to A/" = 8 supergravity amplitudes. The identification of additional cancellations in this the- 
ory [s^ suggests that in four dimensions it may be ultraviolet finite to all loop orders [sS, 
(See also refs. 35|.) String dualities also hint at UV finiteness for A/" = 8 supergravity 
although this is weakened by issues with towers of light non-perturbative states from branes 
wrapped on the co,npact d.„.en=ion. H. Renwkably, computational advances fa gauge 
theory amplitudes can be imported [24] directly into calculations of gravity amplitudes, by 



36|, 



combining the unitarity method with the Kawai-Lewellen-Tye [38[ tree-level relations be- 
tween gauge and gravity theories. This allows cuts of gravity loop amplitudes to be expressed 
as double copies of cuts of corresponding gauge theory amplitudes (2J|. This strategy has 
recently been applied to obtain the three-loop four-point amplitude of A/" = 8 supergravity, 
starting from corresponding A/" = 4 MSYM amplitudes js^ . That computation shows that 
at least through three loops, MSYM and A/" = 8 supergravity share the same ultraviolet 
power-counting. Because they share the same critical dimension for ultraviolet finiteness, 
both are ultraviolet finite in four dimensions. The five-loop planar super- Yang-Mill ampli- 
tudes obtained here will be an important input for obtaining the corresponding five-loop 
supergravity amplitudes. (The supergravity calculation also requires the non-planar contri- 
butions.) 

The paper is organized as follows. In section [Tll we briefly summarize known properties 
of the amplitudes and define the notation used in the remainder of the paper. In section UTTl 



we review the observations of refs. 



3, 



29|, on the exclusive appearance of pseudo-conformal 



integrals in planar four-point amplitudes. Because candidate integrals proliferate as the 
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number of loops increases, we give a systematic procedure in section [IV] for constructing 
these integrals. The results of this procedure at five loops are given in section |Vl along with 
the coefficients of the integrals determined via unitarity. Our ansatz for the amplitude is 
also presented in this section. We then briefly review generalized unitarity in section IVI[ 
A description of the cuts used to determine the integral coefficients is given in sections I VII I 
and IVIIH along with a description of the method of maximal cuts introduced in this paper. 
Our conclusions and some comments on the outlook are given in section IIXI 



II. NOTATION AND REVIEW OF MSYM AMPLITUDES 
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n-2 



We use a standard color decomposition 20|, l39|] for the MSYM amplitudes in order to 
disentangle color from kinematics. In this paper we focus on the leading-color planar con- 
tributions, which have a color structure similar to those of tree amplitudes, up to overall 
factors of the number of colors, N^.. The color-decomposed form of planar contributions to 
the L-loop SU (Nc) gauge-theory n-point amplitudes is, 

^'"Xfi E Tr(TVi,TV.) . . .T'^.(n))A(^)(p(l),p(2), . . . ,p(n)) , (2.1) 
[Any J y 

where Ai^^ is an L-loop color-ordered partial amplitude. We have followed the normalization 

n 

conventions of ref. Here 7 is Euler's constant, and the sum runs over non-cyclic permu- 
tations, p, of the external legs. In this expression we have suppressed labels of momenta and 
helicities, leaving only the indices identifying the external legs. Our convention is that all 
legs are outgoing. This decomposition holds for all particles in the gauge super-multiplet. 
We also define a loop amplitude normalized by the tree amplitude, 



tree 



(2.2) 



Supersymmetry Ward identities [40| guarantee that, after dividing out by the tree ampli 



tudes, MHV amplitudes are identical for any helicity configuration [2l|]. (The complete 
set of these tree amplitudes are tabulated in Appendix E of ref. 2J].) Because four-point 
amplitudes are always maximally helicity violating, this holds for all four-point amplitudes. 
Because it is independent of the position of the two negative helicity legs, M^^^ has com- 
plete cyclic and refiection symmetry. A practical consequence of this is that once a coefficient 
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of a given integral is determined, the coefficient of integrals related by cyclic or reflection 
symmetry follow trivially. 

In our evaluations of four-dimensional unitarity cuts, we use the spinor helicity formal- 



ism 



39 



41| . in which the amplitudes are expressed in terms of spinor inner products, 

--u4ka)[h + k]uAkd) (2.3) 



kb + 



d' 



where u±{k) is a massless Weyl spinor with momentum k and positive or negative chirality. 
Our conventions follow the QCD literature, with [ij] = sign{k'^k^) {j i)* for real momenta 
so that, 

{ij)[ji] = 2ki- kj = Sij . (2.4) 

We also define, 

Xk, = u+{ki), Afe, = u-{ki) . (2.5) 

For complex momenta these two spinors are independent, though they are dependent for 
real momenta. 

In ref. jsl, a conjecture was presented that MSYM amplitudes possess an iterative struc- 
ture, based on an observed iteration of two-loop four-point amplitudes. In ref. j4|, this was 
fleshed out for MHV amplitudes into an explicit exponentiation ansatz to all loop orders. 
Through five loops, the expansion of the exponential gives the iteration relations, 

^1 M^^^ffV^ 



Mf(e) 
Mf(e) 
Mr(e) 

Mf (e) 



kM«(e))' + /(^He) M«(2e) + C(^) + 0(6) 



(e)l ' + M(^) (e) M(2) (e) + (e) M^^^ (3 e) + C^'^ + 0{e) 



(2.6) 
(2.7) 



M«(. 



M«(6)] (e) + M«(e)Mf (e) + \[M^:\ef 



+ 



+ /(^)(e)M«(4e) + C(4) + 0(e), 

M«(^)l (e) - fM«(6)l Mf (e) - M«(e)[Mf ( 



+ M«(e)MW(e) + Mf (e)Mf (e) + f''\e) M«(5e) + C^^) + 0{e 
where f^^\e) is a three term series in e. 



(2.8) 



(2.9) 



.2 AL) 



2 ) 



(2.10) 



and fj;^^ are numbers independent of the kinematics and of the number of external legs n. 
Similarly, the C^^^ are also pure numbers. The constant /d^^ is proportional to the L-loop 



cusp (soft) anomalous dimension 7)^ 



Jo 



1 (L) 



After subtracting the known infrared singularities 
simple exponential form, 



(2.11) 



421 the iteration relation takes on a rather 



J'„(0) = exp 



-7xF«(0) + C 



(2.12) 



where F^^)(0) is the n-point one-loop finite remainder, •jx is the complete cusp anomalous 
dimension, and C depends on the coupling but not on the external momenta. Very recently 
Alday and Maldacena have matched this expression at strong coupling for n = 4 using string 
theory 0. 

The iteration conjecture has so far been confirmed for two- and three-loop four-point 
amplitudes [s], 4| as well for two-loop five-point amplitudes For the four- loop four-point 
amplitude, the integrand is known |l^ and has been shown to generate the correct form 
of the infrared singularities though C(l/e^). This has been used to extract the four-loop 



11| numerically. 



contribution to the cusp anomalous dimension [K 

As the amplitudes are infrared divergent, we need to regulate them. In order to preserve 
the supersymmetry we use the four-dimensional helicity (FDH) scheme 431], which is a 
relative of Siegel's dimensional reduction scheme [4^. 



III. PSEUDO-CONFORMAL INTEGRALS 



Conformal properties offer a simple way to identify integrals that can appear in planar 
MSYM amplitudes [2^ . This observation allows us to easily identify a basis of integrals, 
whose coefficients can be determined via the unitarity method. This greatly simplifies the cut 
analysis because we need to determine only these coefficients to obtain the planar amplitudes. 

Although the underlying theory is conformally invariant, there is as yet no proof that 
only integrals dictated by conformal invariance can appear. One obvious complication to 
providing such a proof is the infrared divergence of the amplitudes, and the subsequent 
need to regulate the integrals (via dimensional regularization), which breaks the conformal 
invariance. As mentioned in the introduction, we therefore call the integrals correspond- 
ing to conformally- invariant ones "pseudo-conformal" . We shall describe in this section 
how to implement this correspondence. In principle, it is possible that individual integrals 
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appearing in the amplitude would have no special conformal properties, yet the complete 
amplitude would retain simple conformal properties because of cancellations between in- 
tegrals. Through four loops 1^, however, only pseudo-conformal integrals appear in the 
four-point amplitude. This provides compelling evidence that this is a general property of 
MSYM four-point amplitudes.^ 

We therefore assume that only pseudo-conformal integrals appear in the five-loop planar 
amplitudes. One way of proving the correctness of this assumption would be to compute a 
sufficient number of cuts in D dimensions to determine the amplitude completely. In this 
paper, we discuss only a partial confirmation. 

Dimensional regulation of the infrared singularities breaks the conformal symmetry. For 
the purposes of exposing the conformal symmetry, we instead regulate the infrared diver- 
gences by taking the external momenta ki off shell and letting the dimension be D = 4. We 
will obtain a pseudo-conformal integral from a suitable conformal integral by reversing this 
change of regulator. 

The authors of ref. j^ol provide a simple way of making manifest the conformal properties 
of planar integrals via "dual diagrams". The dual diagrams provide a direct method of 
identifying all conformally invariant loop integrals.^ In general, conformal properties are 
not obvious in the momentum-space representation of loop integrals, but with a simple 
change of variables encoded by the dual diagrams we can make these properties manifest. 

Let us give an illustrative example. Consider the two- loop double-box integral of fig.[T](a), 

p^ij) — kiY{p — /ci — k2Yq^{q — k^Y^q — /cs — k/^yij) + qY ' 

(3.1) 



where s = {ki + k2Y and t = {k2 + k^y. 

After replacing the regulator as mentioned above, the conformal symmetry can then be 
exposed via the change of variables, 

ki = , /C2 = Xi2, k-i = X23, ki = X34, p = X45, q = Xq4, (3.2) 

where The new variables automatically satisfy momentum conservation 

xai + X12 + X2'i + X34 = <J=^ ki + k2 + k'i + ki = Q . (3.3) 



^ For higher-point amphtudes, the situation is more comphcated, as can be confirmed by checking the 
conformal properties of known results [9] for the five-point amplitudes at tw o loops. 
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The dual diagrams are related to the dual graphs used in graph theory 45|. 
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FIG. 1: The two-loop planar double-box integral (a) and its dual (b) overlaying a faded version 
of (a). In (b) the dashed lines represent a numerator factor of (^24)^x^3 = s^t. This inserted 
numerator factor is needed for conformal invariance of the integral. 



After substituting the new variables into eq. (13.11) and taking D = A, the double-box integral 
takes on a very symmetric form, 



-2\2 4 

I XoA X 



24-^13 



2 2 2 2 2 2 2 

<-y»-" ry-^ <-y»-" ry-^ <-y»^ <-y»^ 



(3.4) 



[.45x15^25x45^35X62^-56 

The conformal-invariance properties follow from examining its behavior under inversion, 



X. 



2 2 ' 

2 J 



(3.5) 



Under this inversion the double-box (13.41) is invariant because each external point 
a^i, a^2, 2:3, X4 appears equally many times in the numerator as in the denominator, while 
the internal points X5,X6 appears exactly four times in the denominator, precisely canceling 
the behavior of the integration measure. The x variables are useful because inversion re- 
spects momentum conservation, which is not true for an inversion of the original momentum 
variables. 

More generally, following refs. [2^, we keep track of conformal weights using dual 
diagrams. To obtain the dual representation, start with the momentum representation 
shown in fig. [It^a) and place internal points 0:5, x^ inside each loop, as well as external points 
Xi, X2, X3, X4 between each pair of external momenta, as shown in fig. [T](b). Following ref. 29 1 
we mark the internal integration points by solid dots at the center of each loop but in most 
cases leave the external points unmarked. Solid lines represent an inverse power of x?^, 
corresponding to the dual propagator l/x^^, which crosses exactly one Feynman propagator 
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h 



h 



FIG. 2: The rung-rule for generating higher-loop integrands from lower- loop ones. 

whose momentum is equal to Xij. Dashed lines represent a positive power of an xjp such 
as the numerator factors of s = 0:24 and t = x^g. The xjj represented by a dashed hne 
correspond the sum of the momenta of the ordinary propagator lines they cross. (The 
dashed lines can be deformed to cross different propagators, but momentum conservation 
ensures that this does not affect the value of the dual invariants x]j.) The dual diagram 
constructed in the way in fig. [H^b), is in direct correspondence to the dual integral (13. 4p . 

We can further restrict the possible set of conformal integrals by requiring that they 
have only logarithmic behavior in the on-shell limit. That is, we are not interested in 
conformal integrals which vanish or diverge with a power-law behavior in any kf, because 
these do not correspond to massless integrals in dimensional regularization. For example, 
numerator factors such as — ^2 allowed. Similarly, factors such as l/x^2 = 1/^2 

are excluded because their power singularities are too severe for the required logarithmic 
behavior of infrared singularities. We then obtain a pseudo-conformal integral, as discussed 
earlier, by replacing the off-shell regulator with the usual dimensional one. 

It is straightforward to generalize this graphical mapping to any loop order, allowing for a 
relatively simple bookkeeping of the change of variables between the momenta and the dual 
Xi variables. The map in eq. (13. 2p exemplifies the convention we use for external momenta 
for all diagrams in this paper. 

The conformal weights are easy to read off directly from the dual diagrams. For a dual 
diagram to be conformally invariant it must satisfy the following: The number of solid lines 
minus the number of dashed lines entering a point Xi must be zero for external points, 
and four for internal points. The conformal weight of four for internal points cancels the 



conformal weight of the integration measure given in eq. (13.51) . As observed in ref. [lOj, a 
consequence of requiring integrals to be conformal is that integrals with triangle or bubble 
subdiagrams are not allowed. 

One class of pseudo-conformal integrals may be understood in terms of the "rung-rule" of 
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FIG. 3: The rung rule maintains conformal weight. If the dual diagram prior to applying the rung 
rule has the proper conformal weight so will the resulting diagram. 



ref. 



231 . This rule instructs its user to generate contributions to an {L + l)-loop amplitude 



from a known L-loop amplitude by inserting a new leg between each possible pair of internal 
legs, as shown in fig.[2J From this set, all diagrams with either triangle or bubble subdiagrams 
are removed. The new loop momentum is integrated over, after including an additional 
factor of (/i + 12)"^ in the numerator, where h and I2 are the momenta flowing through the 
indicated lines. (With the conventions used here it is convenient to remove a factor of i 

n 

from the numerator factor, compared to ref. [23^.) Each distinct contribution should be 
counted only once, even if it can be generated in multiple ways. Contributions arising from 
identical diagrams (that is, having identical propagators) but with distinct numerators count 
as distinct contributions. The diagrams obtained by iterating this procedure are sometimes 
called Mondrian diagrams, because of their similarity to Mondrian's art. 

The rung rule may be understood as a consequence of the conformal properties of the 
integrals. As illustrated in fig. [3], if the starting integral is pseudo-conformal, inserting a rung 
splits the inner loop into two side-by-side loops, and the conformal weight of the central dots 
in the figure is unchanged. However, the upper and lower loop need an additional dashed 
line connecting their central dots to maintain their conformal weight. This dashed line 
corresponds exactly to the factor of (/i + I2Y required by the rung rule in fig. [2j 

The rung-rule, unfortunately, does not generate the complete set of planar integrals 10|. 
However, at least through four loops any diagram that it generates is obtained with the 
correct sign. Here we confirm this observation at five loops, using the unitarity method. To 
obtain the remaining, non-rung-rule contributions to the five-loop amplitude, we start with 
the other pseudo-conformal integrals and determine their coefficients using the unitarity 
method. 



12 




FIG. 4: The four-loop "window" diagram. The hghter colored line running through the diagrams 
is a five-particle cut which separates the diagram into a product of tree diagrams. 

IV. GENERATING THE PLANAR PSEUDO-CONFORMAL DIAGRAMS 

The proliferation of candidate pseudo-conformal integrals with increasing number of loops 
encourages the development of a systematic construction procedure. One approach is sug- 
gested by examining the {L + l)-particle cuts of an L-loop amplitude. Using such a cut we 
can decompose the loop integrals into products of tree diagrams which then simplifies the 
bookkeeping.^ Our procedure will be: 

1. Construct the set of all possible amputated tree configurations on each side of the cut. 

2. Identify all possible loop integrals by sewing each configuration from the left side of 
the cut with each configuration on the right side of the cut. 

3. Identify all possible overall factors in each integral which make it conformal. 

Because the conformal properties are most obvious in the dual representation described in 
section llllt we have found it convenient to work with dual coordinates and translate back 
to the momentum representation at the end. 

A. Constructing all dual tree diagrams 

It is useful to consider first the reverse procedure of splitting an L-loop integral into two 
trees via an l)-particle cut. As an illustration, consider the four-loop "window" integral 

^ At six loops and beyond it turns out that there are pseudo-conformal integrals which do not have an 
{L + l)-particle cut. These can however be obtained from a "parent diagram" containing an (L + 1)- 
particlc cut, by canceling one of the propagators. 
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FIG. 5: The dual diagram corresponding to the window diagram. The hnes with the arrows 
indicate the cut which separates the dual diagram into a product of tree dual diagrams. 

as cut in fig. HI This integral contributes to the four-loop four-point amplitude jlO]. As shown 
in the figure, a five-particle cut separates the integral into a product of two tree diagrams. In 
dual space, we also split integrals into tree diagrams with a five-particle cut. Figure [5] depicts 
the separation of the corresponding dual diagram into tree diagrams, which we carry out in 
two steps. In the first step of the figure we divide the dual diagram along the cut marked 
with arrows, keeping the cut line on both sides. In the next step we drop the lines with 
arrows, giving us dual diagrams corresponding to amputated {i.e. with external propagators 
removed) tree-level momentum-space diagrams. In this representation each dual line crosses 
an internal propagator of the momentum-space diagrams. The diagrams will always have 
a fixed cyclic ordering. The dual-diagram points also respect this cyclic ordering. In this 
example, for the tree amplitude on the left, the points are ordered {x^, Xi, X2, Xg, Xj, Xq, x^}, 
while the points for the tree diagram on the right are ordered {x2, X3, X4, X5, xq, xj, xq}. 

Our systematic construction of all dual loop diagrams simply reverses the pro- 
cess in this example. At L loops we start with two ordered lists of L + 3 points: 
{xi,X2,xl+4,xl+3, ■ ■ ■ ,x^,Xi} foT the left tree and {x2, ^3, X4, X5, . . . , for the right 
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tree. These lists correspond to {L + l)-particle cuts in the Si2 channel; the corresponding 
construction in the S23 channel is easily obtained by relabeling the final result. The as- 
signment of labels xi,X2,X3,X4 to points is determined by the external momenta, and the 
remainder follow from the cyclic ordering. 

We obtain all possible pairs of dual tree diagrams by connecting non-adjacent points with 
1 / xfj propagators in all possible ways such that the lines do not cross. Dual diagrams where 
nearest-neighbor points are connected are not included as they correspond to momentum- 
space diagrams whose external propagators have not been truncated. 

After identifying the possible dual tree diagrams we restore the dual lines representing 
the cut, by retracing our steps in the example shown in fig. [51 That is, in both sets of 
tree diagrams we draw lines with arrows connecting X4 to X5, etc., to Xl+a, which is then 
connected to X2- Once this is done the pairs of tree diagrams can be glued together along the 
lines with arrows, which then gives us the L-loop dual diagrams that we wish to construct. 
At this stage we can remove diagrams trivially related by cyclic or fiip symmetry. 



B. Finding the pseudo-conformal integrals 

Once we have a set of candidate loop-level dual diagrams, we must find the numerator 
factors necessary to make the corresponding integrals conformal. This can be accomplished 



as follows 



10|,l29|: 



If one of the internal points {x^, Xq, . . . , xl^^} appears in less than four dual propaga- 
tors, discard the diagram as it cannot be made conformal. 

To determine possible numerator factors one first identifies all external points from the 
set {xi, X2, X3, X4} appearing in one or more dual propagators, and all internal points 
in the set {xs, Xg, 0:^+4} appearing in five or more dual propagators. All such 
points require numerator factors x'^j to cancel the extra conformal weight. That is, the 
number of times a given external Xi appears in the dual propagators minus the number 
of times it appears in the numerator should be zero. Similarly, for an internal point Xi, 
the number of times it appears in the dual propagators minus the number of times it 
appears in the numerator should be four. To find the conformally invariant integrals 
we sweep through products of all candidate numerators x?- to identify the ones where 
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the conformal invariance constraints are satisfied. (In principle, there might also be 
an overall resulting factor of 1/s = l/x\^ or 1/t = l/a^ia, tiut this does not occur at 
five loops, nor do we expect such contributions to enter the amplitudes with non-zero 
coefficients at any loop order.) 

• If the previous step yields a previously- identified pseudo-conformal integral, go on to 
the next case. Such repeated integrals can arise when a numerator factor cancels a 
propagator or when diagrams are related by symmetries. 

Once we have the set of conformal dual diagrams we can convert these back to momentum 
space with a change of variables, 

ki = X41 , ^2 = 3^12 , ^3 = ^23 , ^4 = 3:34 , 

h=X45, l2 = X56, ... , = a;(L+4)2 , (4.1) 

where the h are the momenta of the lines in {L + l)-particle cut used in the construction. 
Since our construction was only a bookkeeping device for finding pseudo-conformal integrals, 
at the end there is no on-shell restriction on the /j. 

In the next section we apply this procedure to construct a basis of all pseudo-conformal 
integrals appearing in the five-loop planar MSYM amplitudes. 

V. THE FIVE-LOOP PLANAR PSEUDO-CONFORMAL INTEGRALS 
A. The five-loop pseudo-conformal integral basis 

Following the procedure described in the previous section we find a total of 59 inde- 
pendent pseudo-conformal integrals potentially present in the five-loop four-point planar 
amplitude (not counting those related by permutations of external legs). They are shown 
in figs. [71 [HI [9l and [TOl The 'parent' integrals, shown in fig. [TJ have only cubic vertices. The 
remaining integrals have both cubic and quartic vertices. They may be obtained by omitting 
propagators and modifying numerator factors present in the parent integrals. As we shall 
show in the following section using unitarity cuts, the integrals in figs. [Hand [8| appear in the 
amplitude (12.21) with relative coefficients of ±1, which we have absorbed into the definitions 
of the numerator factors in the figures. The remaining ones shown in figs. [9] and [10] do not 
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FIG. 6: The notation used in this section for hsting out the pseudo-conformal integrals contributing 
at five loops. The momentum flow through a line connecting points 5 and 8 gives the momentum 
invariant x^g. 

appear at all. We do not have an explanation for the remarkable simplicity of the coefficients 
of the integrals, but presumably it is tied to the superconformal invariance of the theory. 

We draw the diagrams in momentum space, but also include the relevant Xj for tracking 
numerator factors. The numerators are written out as Mandelstam variables s = {ki + ^2)^ 
and t = {k2 + ^3)^ or as dual invariants, x'^j . As discussed in section [XTTl a dual invariant 
xjj is equal to where K is the total momentum flowing through a line spanned between 
points i and j. For example, in fig. [6|, = (/i — h + h)^ = xf^- 

Some of the integrals have identical sets of propagators but differing numerator factors. 
These sibling integrals would be identical were one to omit the numerators. Examples are In 
and /12 or I21 and in fig. [71 The numerator factors often have different symmetries than 
the propagators in any given integral. The different numerator factors in sibling integrals will 
also typically have different symmetries. For example, integral /22 is completely symmetric 
under a cyclic permutation of its arguments, {1, 2, 3, 4} {2, 3, 4, 1} (corresponding to a 7r/2 
rotation of the diagram in fig. [7]) and under fiips, {1, 2, 3, 4} ^ {4, 3, 2, 1} (corresponding to 
refiection of the diagram). Its sibling /21, in contrast, has only one symmetry, {1, 2, 3, 4} — »• 
{3,4,1,2} (corresponding to a rotation of the diagram by vr radians). Accordingly, J22 
appears only once in the amplitude, but /21 appears four times. This makes it inconvenient 
to combine them into a single integral. 

All planar five-loop pseudo-conformal integrals with the exception of I55 have at most 
four-point vertices. (J55 has quintic vertices where the external legs attach. Internal quintic 
vertices do not occur in conformal integrals until seven loops.) 



17 



7. 



(/l) 

2. 







7. 







5t X45X3g 





6. 








9. 










5. 


6. 





01 -^45 



,y,2 ^2 



5. 


6. 








9. 



ih) 



q2j-™2 ™2 ™2 
* ""^as-^ 17-^38 



(^7) 



T"^ -7^2 ™2 

X35X28J'47 



ih) 



ih) 



{I 



4. 



C^^T*^ /y«2 ^2 

56 J^27-^37'^46 





6. 






7. 









4. 

(/13) 















6. 




5 

• 


7. 


3 

• 


q2j-™4 ^2 

* tXggXgg 


5 

• 




















9. 



4. 

(As) 



524^2 ,yi2 ™2 
t.X3gX4gX29 





6 

• 




5 

• 




9. 



4. 

(/17) 
2 



it X59 



qV^3q^2 ^2 
it X4gX37 



(/lo) 



S2j^^2 ^4 





6 

• 


7 

• 





(/14) 



Q^-f-zy'^ ^2 

i tX59X3g 





6 






• 




5 




9 


• 




• 



ihs) 



C^/'T*^ ry«2 ™2 

i 6X35X45X37 





6. 






5 






• 


7 






• 





4. 

(/19) 



^» (-j24-™2 ™2 ^2 
i t.X35X4gX5g 





6 






• 


5 






• 




9 






• 



^2 ™2 ™2 
it X28X45X5g 



5 


1 6 


• 


X * 




y9 


• 


1 • 



(/20) 



4. 



524^2^2 ^2 




FIG. 7: All pseudo-conformal integrals with only cubic vertices that contribute to the amplitude. 
The relative signs are determined from unitarity cuts in sections IVIII and IVIIII 

B. The five-loop four-point amplitude 



In sections I VIII and IVIIIt we evaluate a sufficient number of cuts in order to determine the 
numerical prefactors of each pseudo-conformal integral as it appears in the amplitude (12.21) . 
We find that the complete five-loop four-point MSYM planar amplitude is, 
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FIG. 8: All pseudo-conformal integrals with cubic and quartic vertices that contribute to the 
amplitude. The relative signs are determined from unitarity cuts in sections IVIII and IVIIIl 



+ 2/n + 4/12 + 4Ji3 + 4/14 
+ 2J21 + 2/23 + 4/24 + 4/25 



4/15 + 2/16 + 4/17 + 4/18 + 4Ji9 + 4/20 
4/26 + 2/27 + 4/28 + 4/29 + 4/30 



2/31 + + 4/33 + 2/34 + {s^ t}) 



(5.1) 



where the integrals are shown in figs. [7] and [H and Mf^ is defined in eq. (12. 2p . There 
are a total of 193 integrals in the sum. As the integrals depend only on the kinematic 
invariants s and t, instead of having leg labels, each integral can appear only as Ij{s,t) or 
as Ij{t,s). In eq. (15. ip . we have suppressed the arguments "(s,t)" and combined identical 
terms, leaving a symmetry factor in front. The relative signs between integrals, determined 
from the unitarity cuts in sections IVIII and IVIIIt have been incorporated in the numerator 
factors in figs. [7| and [H though we have chosen to leave an overall sign outside the integrals. 
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FIG. 9: A class of pseudo-conformal integrals which do not contribute to the amplitude, as deter- 
mined from the unitarity cuts. All these have exactly one factor of st. 

The normalization factor of 1/32 follows the conventions of refs. and accounts for the 

factor of 2^ in eq. (12.11) . The integrals in eq. (15.11) are therefore normalized as: 

5 



(5.2) 



where the li are five independent loop momenta, A^ is the numerator factor appearing as the 
coefficient of the diagrams given in figs. [TlfTOl and the correspond to the propagators of 
the diagrams. 

To understand the relative signs of the diagrams we classify terms into those derived 
from the rung rule and the rest. Any diagram generated by the rung rule in fig. [2] simply 
inherits the sign of the lower loop diagram from which it was derived. This gives the 
correct numerator factor, including the sign, for all contributing integrals containing only 
cubic vertices, except for I21 - the only diagram in fig. [7] having a non-rung rule numerator. 
Integrals with quartic vertices, such as J24, Irj-, -^28 and 129, are given by the rung rule apphed 
to known |lO| four- loop diagrams. Using two-particle cuts, other examples of integrals whose 
prefactors are easy to understand are /23, /25, 12%-, hi and 142- The latter two have vanishing 
coefficient because their four-loop parent diagrams also have vanishing coefficients. Integral 
/33 can be understood in terms of a rung inserted between an external leg and an internal 
line. 

Another class of prefactors and signs can be understood from a "substitution rule". 
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FIG. 10: The non-st class of pseudo-conformal integrals. They do not contribute to the amplitude. 




FIG. 11: The "substitution rule" is a rule for replacing any four-point vertex with terms in a four- 
point amplitude. Here integral I22 is obtained by substituting the pseudo-conformal box integral 
appearing in the one-loop amplitude into the central vertex of the four-loop "window" diagram. 



Consider diagram I22 in fig- 13 As shown in fig. [HI this diagram inherits its prefactor and 



21 



sign by replacing the four-point vertex by a one-loop box integral. (The one-loop box enters 
with a relative plus sign.) The numerator factor XggXgg of the substituted box is simply the 
factor needed to make the box conformal. The negative sign is inherited from the sign of 
the four-loop diagram. Also other signs can be understood from this substitution rule. For 
example, the sign on /30 and the zero coefficient on 153 follow from similar substitutions on 
four loop conformal diagrams. This rule can more generally be understood as a substitution 
of the normalized four-point function into a four-point vertex, which can obtained using 
generalized cuts. 

With the rung rule, two-particle cuts, and substitution rule we may understand the signs 
of all diagrams that appear in the amplitude, except for /31, 1^2 and J34. As of yet we have 
not found a rule giving the sign of these diagrams, other than resorting to computations of 
cuts. 

At four loops [1^, integrals not containing at least one factor of s and also one factor of 
t are absent from the amplitude. This is again true at five loops. Factorization arguments 
using complex momenta can give a suggestive explanation of this property. As we already 
noted, supersymmetry identities ensure that after dividing by the tree amplitude, the non- 
vanishing MSYM four-point loop amplitudes are identical for all external helicity and particle 
configurations. Consider then the helicity configuration 1~, 2"^, 3~, 4+, with all external legs 
gluons. The tree amplitude 

^-'(12)(23)(34)(41) - ' St ' ^^-^^ 

factorizes in both the s and t channels into products of three-point vertices. The absence of 
compensating factors of s and t would imply factorization into one-particle irreducible loop 
three vertices. Such vertices have not appeared in the factorization of any previous MSYM 
amplitude, and so it is not surprising that the offending integrals, shown in fig. [10], do not 
contribute here either. 

We may understand the remaining vanishing coefficients using the known harmonic- 
superspace power counting of MSYM 46||. It is compatible in D = 4 with six powers 
of loop momenta canceling from integral numerators. The missing engineering dimensions 
of the amplitude are then supplied by external momenta, requi ring at least three powers of 



either s or t. This result agrees with the arguments of refs. 
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24l |. which provide a bound 
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on dimensions for which the L-foop amphtude is ultraviolet finite, 

D<1 + ^, iL>l). (5.4) 

It is natural to assume that this power counting holds independently for each integral. This 
rules out integrals which do not have at least three powers of s or t. This includes all of 
those in fig. [9] and those with either a single power of or in fig. [TOl 

C. All-loop structure 

Inspecting the contributions of the integrals in the basis to the five-loop four-point am- 
plitude reveals the following features: 

• All pseudo-conformal integrals containing a factor of s'^t or t'^s (and possibly additional 
powers of s or t) enter the amplitude with relative weight of +1 or —1. 

• Any pseudo-conformal integral without a factor of s^t or t'^s has a vanishing coefficient. 

• All integrals that could be obtained from the rung rule or from two-particle cuts inherit 
their weights from the lower-loop integrals used to construct them. 

• All contributing integrals satisfy the ultraviolet finiteness bound (15. 4p . 

These observations seems to suggest that, in general, the set of pseudo-conformal integrals 
with nonvanishing coefficients are the ones with a prefactor divisible by either s'^t or t^s. 
However, at six loops a new structure appears where a pseudo-conformal integral is a simple 
product of lower-loop integrals, as displayed in fig. [121 We have checked that the conformal 
integral in fig. [12] does not contribute to the amplitude, although its prefactor is divisible by 
sH. While this particular integral does not appear in the amplitude, its existence suggests 
that at higher loops there will be additional classes of pseudo-conformal integrals with 
vanishing coefficients. 

VI. GENERALIZED UNITARITY 



The unitarity method 17], [l8|, [l9|, |20|, |2l|, [22, |26| has proven an effective means for com- 
puting scattering amplitudes in gauge and gravity theories. So-called generalized unitarity 
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FIG. 12: A pseudo-conformal integral with a vanishing coefficient in the six- loop amplitude. 



is particularly powerful for computing amplitudes 22, 1261], as it allows an L-loop ampli 



tude to be bui 
momenta 26, 



t directly from products of tree amplitudes. When combined with complex 



27 



47( 1 . it allows the use of maximal cuts, in which all propagators in an 



integral are cut. (The term "generalized unitarity", corresponding to leading discontinuities 
of diagrams, dates back to ref. j48|.) 

We begin our discussion with a brief review, including earlier applications of maximal 
cuts to the computation of two-loop amplitudes 26|, l28|]. We record a number of observations 
useful for computation at higher loops. In section lVI Al we modify the maximal-cut procedure 
and use it to determine the coefficients of all pseudo-conformal integrals appearing in the 
MSYM five-loop four-point amplitudes efficiently and systematically. 



Maximal cuts 



Cut calculations can be simplified by increasing the number of cut legs. This isolates 
a smaller number of integrals, making it simpler to determine the values of their coeffi- 
cients. This technique is especially powerful for computing one-loop MSYM amplitudes, 
because only box integrals can appear [l7i]. As observed by Britto, Cachazo and Feng 26|. 
taking a quadruple cut, where all four propagators in a box integral are cut, freezes the 
four-dimensional loop integration. This allows its kinematic coefficient to be determined 
algebraically , with no integration (or integral reduction) required. The use of complex mo- 
menta, as suggested by twistor space theories [271], makes it possible to define massless three 
vertices and thereby to use quadruple cuts to determine the coefficients of all box integrals 
including those with massless external legs. 

For three massless momenta ka, kf, and k^ = —{ka + k^) one has the following consistency 
requirement 



= k^c = (ka + = 2ka ■ k, = {ah) [6. 



(6.1) 
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For real momenta in Minkowski signature, A^^ and A^^ (see eq. (12 .Sp ) are complex conjugates 
of each other (up to a sign determined by incoming or outgoing nature of the corresponding 
particle). Hence if (ab) vanishes then [ab] must also vanish. This constraint holds for all 
three legs a, b, c, leaving no non- vanishing quantities out of which to build a three vertex. If 
the momenta are taken to be complex, however, the two spinors A^^ and A^^ are independent. 
This gives two independent solutions to eq. (16. ip . 

{ab) = 0, or [ab] = 0, (6.2) 

with the other spinor product non-vanishing in each case. In the three-gluon case, there 
are overall two possible solutions: all As proportional, and hence all (ij) vanishing, or all A 
proportional, and hence all [ij] vanishing. This means that exactly one of, 

4"^ = Ar>-, h\ c+) = -^^4t-T ' (6-3) 

[ab\ [ca\ 

At^ ^ Ar{a\b-,c~) = I , (6.4) 

{ab) {ca} 

does not vanish. Similar statements hold for amplitudes involving fermions or scalars: one of 
the two independent helicity configurations will not vanish. The non-vanishing amplitudes 
involving a fermion pair are, 

/ir(a7,6;,c+) = -zM!, (6.5) 



Ar%a|,bf,c-) = ~^^, (6.6) 



where the subscript / denotes a fermionic leg. Similarly the non-vanishing scalar amplitudes 
are, 

4^^^(a;,&^c+) = -.MM, (6.7) 

Ar{atA,c-')=^^-'^^^, (6.8) 

where the subscript s denotes a scalar leg. (For complex scalars, the two helicities correspond 
to particle and antiparticle.) We have chosen the signs in these amplitudes to be consistent 



with the supersymmetry Ward identities 



conjugation rules of ref. 49 1 



3, 



as given in ref. [18[], and with the parity 



The method of quadruple cuts has been generalized by Buchbinder and Cachazo 



28] to 



two loops using hepta- and octa-cuts. Although the two-loop four-point double-box integral 
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FIG. 13: A pictorial representation of two kinematic solutions to the four-point hepta-cut equations. 
The diagrams representing the remaining four solutions can be obtained by reflection symmetries 
of these. A '@' vertex represents a three-point tree amplitude involving only A spinors and a 
'0' vertex represents an amplitude involving only A spinors. All lines are cut and carry on-shell 
momenta. 

only has seven propagators it secretly enforces an additional, eighth constraint, yielding an 

] point 



octa-cut which localizes the integration completely. But as the authors of ref. 
out this last cut condition is not really necessary for the evaluation of the two-loop four- 
point amplitude: the integrand was already independent of loop momenta after imposing 
the constraints from the seven on-shell propagators in the hepta-cut. We will return to this 
point below. 

Let us then focus on the constraints imposed by the delta functions corresponding to 
the propagators alone. In the hepta-cut construction of ref. [28|, various classes of solutions 
are allowed by the seven delta function constraints arising from localizing the propagators. 
In the context of generalized unitarity these delta functions correspond to solving the on- 
shell cut conditions If = 0. As discussed above the solution to these conditions is always 
complex when three-point vertices are present. These cut conditions have a discrete set of 
solutions because of the two-fold choice in eq. (16. 2p at each three-point vertex. Each of 
these solutions depends on continuous parameters, corresponding to the degrees of freedom 
not frozen by the cut conditions. The discrete choice coincides with the choice of three- 
point amplitude at each vertex A'^^ as given in eqs. (16.31) and (16. 4p (or similar vertices 
for fermionic and scalar lines), which suggests a convenient way to represent the possible 
solutions using additional labels at the vertices of the cut diagrams. For the four-point hepta- 
cut two inequivalent arrangements of three-point vertices are shown in fig. [131 External legs 
represent outgoing external momenta while internal lines represent cut propagators and thus 
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on-shell loop momenta. The signs inside the blobs in the diagrams indicate the corresponding 
choice for the three-point vertex, and implicitly, that the spinors of the opposite helicity are 
proportional. A '0' vertex will have all A spinors proportional to each other, so that the 
vertex is built out of A spinors of the attached legs, while the roles of the two kinds of 
spinors are interchanged for a '0' vertex. We can sum over all possible solutions to obtain 



the multiply-cut integrand, as was done in ref. [28|]. For the sevenfold cut of the double-box 
diagram, there are six distinct solutions of the two types shown in fig. [121 

Once we choose external helicities, as in fig. [TU the blobs then dictate the possible 
assignments of helicities for internal lines. The rules for finding the complete set of kinematic 
solutions associated with a given assignment of plus and minus labels to a diagram are as 
follows: 

• A '0' label means the three A spinors corresponding to the lines attached to the 
blob are proportional to each other. Similarly, a '0' label denotes having the three A 
spinors proportional to each other. 

• If one of the lines attached to a vertex is an external line ki then the spinors are 
proportional to an external spinor, either A^- or A^-. 

• If a '0' vertex is directly connected to another '0' vertex then all As of the lines 
attached to both vertices are proportional to each other. A similar statement holds 
for the As of two connected '0' vertices. 

• If there is a chain of vertices of the same sign connecting any two external lines then the 
diagram vanishes, because one cannot solve the on-shell and momentum conservation 
constraints for the diagram. A solution would require that two external spinors of the 
same type are proportional to each other; this cannot be true in general, because they 
are independent. 

Applying these rules to the four-point double-box does indeed give the six allowed solu- 
tions of the two types shown in fig. [131 The remaining solutions are related to the depicted 
ones by flip symmetries. The complete set of solutions to the cut constraints is solely 
determined by the topology of a given diagram. Each solution determines a pattern of 
'0' and '0' vertices in the diagram. For each solution to the cut constraints, one of the 
two types of three-point amplitudes vanishes at each vertex. This pattern (along with the 
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FIG. 14: A singlet hepta-cut (a) and one of the two helicity configurations (b) in the simplest 
non-singlet cut. The latter allows gluons, fermions and scalars to propagate in the loop indicated 
by a dashed circle. The other configuration is obtained by flipping all the helicity signs of the legs 
in this loop. All lines in this figure are cut and carry on-shell momenta. The arrows in (a) refer to 
the direction of momentum flow. 

topology) will in general restrict the helicity assignments along internal lines, and may also 
restrict the particle types allowed in different internal lines. 

The strongest constraint that can arise in a kinematic solution is the restriction to a 
single allowed helicity configuration for the internal lines. We will refer to this configuration 
as a "singlet" . In this case only gluons can propagate inside the diagram, as in fig. [T4l(a) . 
Fermions or scalars are not allowed because the only potentially non-vanishing vertices are of 
the wrong type and vanish for the given solution. The second- strongest constraint allows two 
helicity configurations. In such configurations the particle content is purely gluonic except 
for one loop in which any particle type can propagate, as shown in fig. [TW b). (A fermionic 
loop always allows two helicity assignments, corresponding to interchanging fermion and 
antifermion, and the same is true for complex scalar loops.) 

Solving for the spinors in any diagram is then straightforward. Consider the singlet case, 
diagram (a) in fig. [141 The on-shell conditions together with momentum conservation at 
each vertex give a set of equations that must be satisfied. 
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c<: oc Ai3 , A^A/4 — A/gAj^ — A^gA^g , 

A/4 oc Aig oc A/2 , A^A/4 = A/jA/g — A/jA/j 



The solution to these equations is, 



A/. 


= Afci , 




A/2 


= ^Afcj 


-A, 


A/3 


= Afc^ 




K 


= A/2 




A/s 


= A/2 




A/a 


= Afc, ^ 


- A/2 


A/7 


= A/j4 





A/i = ^Afc2 , 

fc2 5 A/2 = Afcj , 

A/a = ^Aa;2 + Afei 



[2 3] 



Au = A.a (6-10) 

A/s = Afca + A/4 

A/„ = A/, 



\ -\ Ml 

^/7 ^fc3 ^/2 4) 

where ^ is an arbitrary parameter, corresponding to the remaining degree of freedom in 
the integration not frozen by the hepta-cut. Since a bispinor p"" = A'^A" is invariant under 
a rescahng of the spinors, (A", A") (/5A" ,/5~^A"), the above solution can be written in 
many other forms. In addition, there is a choice as to where to include the remaining 
degree of freedom. While individual three-point amplitudes ^43^^ are not invariant under 
this transformation, the product of amplitudes forming the cut is invariant. 



B. Solving for integral coefficients using maximal cuts 

We now consider how to solve for the coefficient of an integral using the maximal cuts. 
At two loops, there is only a single conformal integral, the double box. It can appear, of 
course, in both s- and t-channel configurations, but the hepta-cut shown in figs. [13] and [H] 
selects only the s-channel double box. Our candidate expression for the amplitude is then, 

A(2)(1,2,3,4) = cA*'-"^(l,2,3,4)/(')(s,t), (6.11) 

where I^'^\s,t) is the pseudo-conformal two-loop double-box integral in eq. (13.11) and c is a 
coefficient that we need to solve for. This integral contains a factor of s'^t in the numerator, 
which as we shall see is necessary for satisfying the cut conditions. 
Imposing the sevenfold cut condition, we obtain, 

cs2tA*'"^^(l,2,3,4) / d\d%l[6{l^) 

i=i 

= I [ d\d% n m E(^rir^(r^f3r^(4r^(5r^(6r)/. , (6.12) 
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where A^^^*^ is the three-point amphtude corresponding to one of the six three vertices and 
the sum over hehcities h runs over all possible helicity and particle configurations. (We 
have taken the loop integrals to be four-dimensional for the purposes of our discussion here, 
but in any explicit evaluation of the integrals they should be continued to D dimensions to 
regulate the infrared singularities.) 

As discussed above, the delta-function constraints are solved by a discrete set of solutions, 
so we obtain 

csHAT%1,2,3,4) J d\d% j d^j2Jj^\^i-^f')S\k-lf') (6.13) 



; J d%d% I de E - if') ^% - if') E (Afif ^sr^faf ^Sf 4r4 



1 tree\ 
^(6) )j,h , 



where j runs over the different kinematic solutions, f"^^ and 1^°^^ are the values of the 
independent loop momenta expressed in terms of the external momenta, and the remaining 
degree of freedom is ^. For each discrete solution j, only a subset of helicity and particle 
configurations denoted by Hj gives a non-vanishing contribution. The Jacobian from the 
change of variables is Jj. 
Buchbinder and Cachazo 



28| noted that the integrand is constant after imposing the seven 
cut conditions arising directly from cutting propagators, without need to impose the eighth 
cut condition. Another curiosity they noted is that all six discrete kinematic solutions for 
the hepta-cut give the same answer for the amplitude. This was true for kinematic solutions 
that permitted only gluons in the two loops as well as for those which also permitted fermions 
and scalars. This simplicity is related to the absence of terms which integrate to zero upon 
performing the loop integral.^ 

We will exploit this observation, and assume its generalization to higher loops. It allows 
us to match integrands, and indeed to pick individual solutions to match the left-hand and 
right-hand side of eq. (16.131) . determining the overall coefficient. That is, we assume that 
there is a single overall coefficient c to solve for in front of each integral, instead of a different 
contribution for each solution. This also avoids any need for integral reductions or analysis of 
the integrals, and translates integrands into algebraic coefficients of integrals. Our knowledge 



^ This property is special to four-point amplitudes and is already violated at one loop for five-point ampli- 
tudes. 
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of an integral basis — given by the pseudo-conformal integrals — is not essential but greatly 
simplifies the extraction of these coefficients. This equality of contributions from different 
solutions is likely special to MSYM at four points or perhaps to conformal supersymmetric 
gauge theories more generally. In general, there is no reason to expect solutions which allow 
different particle types to circulate to yield equal answers. 

The assumption can be checked directly, of course, by comparing different solutions. 
While we have not checked it exhaustively, it does pass the large number of such compar- 
isons that we have carried out. The use of the assumption and the maximal-cut procedure 
described here also leads to a determination of coefficients at three and four loops in agree- 

nn 

ment with known answers [J, [lOj. Furthermore, a violation would likely lead to inconsistent 
determinations of integral coefficients at five loops; we find no such inconsistency. We can 
also rely on cross checks from non-maximal cuts. 

(The reader may be puzzled by the appearance of complex solutions in what was originally 
an integral over real loop momenta. This is not special to the amplitudes under consideration 
here. In extracting the cut by replacing propagators with delta functions, one must sum over 



complex solutions as well as real ones 



26| . This was necessary in other circumstances such as 



evaluating the connected prescription for tree-level gauge-theory amplitudes 50| in twistor 



string theory [27|. It can also be understood by reinterpreting [51[ the original integral as 
a fourfold contour integral in each component of the loop momentum, and replacing the 
propagators by products of an expression of the form [27ri(/f — times Jacobians; 

Cauchy's theorem makes it act like a delta function, but allowing complex solutions. The 
details are not important to us here because we are only determining coefficients and not 
evaluating any integrals.) 

Choosing one of the kinematic solutions then gives, 

csHA'ni, 2, 3, 4) = z 5:(4f Agf Agf Agf Agf /Ifef ). , (6.14) 

h 

where h runs over the helicity configurations and particle content with non-vanishing con- 
tributions for the given solution. 

Obviously, it is advantageous to choose the simplest solution, where the kinematics re- 
stricts us to the fewest possible particle types circulating in the loop. The best choice is a 
singlet where only gluons contribute. Using the kinematic solution f l6.10p . corresponding to 
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FIG. 15: A singlet iterated two-particle cut of the two-loop four-point amplitude. 



the singlet solution in fig. [TWa). we have, 

/- e^j. 4^66/1 - 9+ q+ /I _ • 4tree /(tree /(tree /(tree /(tree /(tree ( P, ^ ^\ 

CS l/i^ [L ,Z ,6 ,4 ) — ^(2) ^(3) ^(4) ^(5) ^(6) ' 1,0. iO J 

with no sum over intermediate helicities. In the singlet case there is only one term in the 
helicity sum and all six three-point amplitudes are purely gluonic. Here, represents one 
of the three-gluon tree amplitudes in eqs. (16. 3p and (16. 4p . with the plus and minus labels on 
these amplitudes matching the labels of the vertices in the figure. We have confirmed that 
eq. (16.151) holds for any value of the arbitrary parameter ^ (other than ^ = 0, where the 
right-hand-side of eq. (16.151) is ill defined) in the kinematic solution (16.101) . This equation 
then determines c = -|-1, so that the pseudo-conformal double box integral in fig. [T] appears 
in the two-loop amplitude with a coefficient of -|-y44'''^'^(l~, 2+, 3+, 4~), in agreement with 



known results 



23 



24|. 



C. Generalized unitarity with real momenta at two loops 

The kinematics in maximal cuts is highly constrained. It is therefore useful to have a way 
of checking results using less-restricted kinematics. Let us begin by considering generalized 
(non-maximal) cuts which are well defined for real momenta in four dimensions. 

As an example of a four- dimensional generalized cut, consider the two- loop iterated two- 
particle cuts shown in fig. [151 This helicity configuration has the property that it is a singlet 



under supersymmetry transformations [2^ , the only contributions coming from gluon inter- 
nal states. Hence we call it the "singlet" contribution. The remaining contributions, contain- 
ing all other helicity and particle-type assignments, we will collectively call the "non-singlet" 
contributions. These latter contributions transform into each other under supersymmetry. 
(The action of supersymmetry on the A/" = 4 amplitudes is described in appendix E of 



ref. |2^.) 
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FIG. 16: A singlet contribution to the two-loop three-particle cut. Only gluons enter in the loops. 



The singlet configuration is especially simple to evaluate because it involves only a single 
particle type, similar to the singlet configurations of the maximal cuts. At higher loops, the 
number of different particle and helicity configurations grows rapidly. If possible, it would 
be simpler to use only singlet configurations to confirm our ansatz for the amplitude, just 
as with maximal cuts. Unlike maximal cuts, however, the generalized cuts considered here 
do not enforce a particular choice of internal-line helicities or particle types. Nonetheless, 
in special cases, the singlet can be used to determine the coefficient of integrals in the 
amplitude. 

For example, in the iterated two-particle cuts at two loops shown in fig. [151 the singlet 



contribution gives exactly the coefficient of the double box integral [23|]. With this external 
helicity configuration, this cut has no non-singlet contribution. The non-singlet contribution 
appears in the other channel and gives the identical result. In the three-particle cut, however, 
the singlet and non-singlet contributions appear in the same cut and are not identical. This 
cut does nonetheless have simple properties that we can exploit. The singlet contribution. 



depicted in fig. [161 has been previously evaluated in refs. [23|, |2J], with the result 



^singlet ^ ^tree(^-^2-,/3+,/+ /^)xAf^(3+,4+,-/r,-/2",-^3") 

^ r .2 tr+[l/i43/32] 

where tr-|-[l/i ■ ■ ■] = tr[(l ± 75)1^1/1 ■ ■ ■]/2. The tree amplitudes that appear are, 

/itree/1- 9- 7+ 7+ i+\ _ ■ ^) 

Ar{s\A\ -/r, -/,-, -/3-) = -^[3 4] [4 HO] [Hi) H^^H^) Hs)] [H3)3] • ^^-'^^ 

The non-singlet contribution arising from the contribution of all other helicity and particle 
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configurations crossing the cut is a bit more compficated to evaluate, and is equal to [23|, ) 



^non-singlet ^ _ / I fS4l^ tr, [lZi43/32] 

^ M^i + ^2)^(^2 + /3)^(/3-A:3)2(/i-A:4)2(/3 + fc2)2(/i + A:i)2- ^ ■ ^ 

The 75 terms in the singlet and non-singlet appear with opposite signs. In the sum over 
singlet and non-singlet contributions to the cut, the 75 terms therefore cancel algebraically at 
the level of the integrand. (Alternatively, the difference between the singlet and non-singlet 
contributions integrates to zero.) From a practical standpoint, it is easier to compare cuts 
with target ansatze prior to integration and to use only the singlet, so the key observation is 
that we may use is that the non-75 term of the singlet is exactly half the total contribution 
to the cut. 

We need to extend these observations to higher loops in order for them to be useful. We 
have confirmed that the same properties hold at three and four loops for any combination of 
two- and three-particle cuts composed only of four- and five-point tree amplitudes. (If six- or 
higher-point tree amplitudes are present in the cuts, non-MHV configurations with three or 
more negative and three or more positive helicities enter into the computation, which renders 
the structure of supersymmetric cancellations more elaborate and prevents us from using 
the singlet contribution alone to evaluate the cut.) We will assume that this observation 
continues to hold true at five loops. With this assumption we will be able to check (in 
section IVIII|) the coefficients of a variety of pseudo-conformal integrals, using only singlet 



cuts. This is a strong consistency check, because it relies not only on the coefficient under 
examination being correct, but also on the assumption remaining valid. (It seems extremely 
implausible that a breakdown of the assumption at five loops could be compensated by an 
incorrect coefficient.) 

To do better we need to sum over all states crossing the cuts. Moreover, a proper 
treatment of the cuts requires that the cuts be evaluated using D-dimensional states and 



momenta 



19 



2Q|, 



211]. This ensures that no contributions have been dropped, as can hap- 
pen when four-dimensional momenta are used. At one loop, the improved power counting 
of supersymmetric theories allows one to prove a theorem that unitarity cuts with four- 



dimensional momenta are sufficient to determine dimensional 
amplitudes (that is, "near" four dimensions) completely 



17 



ly-regulated supersymmetric 



18|. (The regulator must of 



course maintain manifest supersymmetry; as mentioned earlier, we use the four-dimensional 
helicity scheme (FDH) to do so. In this scheme, the helicity algebra is always four dimen- 
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sional, but the momenta are continued to D = 4 — 2e dimensions.) Unfortunately, no such 
theorem is as yet known beyond one loop. A subtlety in deriving such a theorem arises 
from infrared singularities: the singularities in one loop can effectively "probe" the 0{e) 
contributions from another loop, the product giving a surviving contribution even as e ^ 0. 
When computing with D-dimensional momenta, one can no longer use the spinor helic- 



ity representation 4l|, which makes expressions for tree amplitudes used in the cuts more 
complicated. A good way to ameliorate this additional complexity is to consider instead 
A/" = 1 in ten dimensions super- Yang-Mills dimensionally reduced to D = A — 2e dimensions. 



The remaining states are completely equivalent to those of MSYM in the FDH scheme 
except that the bookkeeping of contributions is much simpler. 



At two loops all cuts of MSYM amplitudes were evaluated in D dimensions 521], providing 



a complete proof of the planar and non-planar expressions for the MSYM amplitudes first 



obtained in ref. 23| using four- dimensional momenta. At three loops, we have also re- 
evaluated the planar amplitude using D-dimensional cuts. The four-loop planar amplitude 
has been evaluated using D-dimensional cuts, assuming that the full result (as an abstract 
tensor in polarization vectors and momenta) is proportional to the tree amplitude, and 
making the reasonable assumption that no contributions can have a triangle subintegral. 
(Only the terms involving polarization vectors dotted into each other, after tensor reductions, 
were evaluated explicitly. Also, one can rule out all bubble and some triangle subintegrals 
using supersymmetry along with generalized unitarity.) In section IVlIIt we will make use 
of these results to provide non-trivial evidence in favor of the various assumptions we have 
used to obtain our ansatz for the five-loop four-point planar amplitude. A complete proof 
would require additional £)-dimensional cuts be evaluated in order to confirm the coefficient 
of every potential integral that might appear, including non-conformal ones. 

VII. MAXIMAL CUT TECHNIQUE FOR DETERMINING INTEGRAL COEFFI- 
CIENTS 

A. Overview of maximal cut method 

In this section, we further develop the maximal-cut technique for higher loops using the 
observations of the previous section. This allows us to extend the two-loop maximal cuts of 
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Buchbinder and Cachazo [28| to higher- loop orders. Unhke Buchbinder and Cachazo, we do 
not require the loop integration be frozen by the cut conditions. That is, we do not require 
the number of cut conditions to match the number of loop integrations. As discussed in the 
previous section, instead, we perform all evaluations of the cuts at the level of the integrand, 
prior to performing any loop integrations. 

Moreover, we do not solve for all possible kinematic configurations satisfying the cuts. 
We instead focus on those solutions which allow the simplest determinations of the coeffi- 
cients of the pseudo-conformal integrals as they appear in the amplitude. As discussed in 
section IVI At the simplest kinematic solutions are the singlets, to which only gluons con- 
tribute. At five loops it turns out that only four pseudo-conformal integrals do not have 
singlet solutions. However, even in these cases one can choose kinematics which forces the 
fermion and scalar contributions into specific loops, again greatly simplifying the determi- 
nation of the coefficients. In order to speed up the process of extracting the coefficients we 
solve the constraint equations numerically, although in some cases we find it useful to solve 
the constraints analytically. 

At two-, three- and four-loops, where the complete results for the four-point planar ampli- 
tudes are known {4, 10, 23 1, we have confirmed that singlet maximal cuts correctly determine 



the coefficients of all integrals appearing in the amplitudes. This suggests that the same will 
be true at five loops. Again, if this were not true it would reveal itself as an inconsistency 
in the results. In particular, we would find that different kinematic solutions of the cut 
conditions would lead to inconsistent determinations of integral coefficients. We would also 
find inconsistencies with cuts with less restrictive kinematics. 

A drawback of our maximal-cut method is its reliance on four-dimensional spinor helicity, 
which as mentioned in section IVI CI might drop contributions. Nevertheless, it does provide 
a relatively simple and systematic means to obtain an ansatz for the coefficients of integrals 
that appear in the amplitude. 

To evaluate the coefficient of any of the five-loop integrals with only three-point vertices 
shown in fig. [7l we cut all 16 propagators. Similarly, the coefficient of integrals with quartic 
vertices in fig. El can be obtained by cutting all of the propagators present. This is of course 
fewer than the number of propagators present in diagrams with only cubic vertices, so some 
of the integrals with only cubic vertices can contribute to the cut. We must therefore subtract 
out all such contributions, to obtain the coefficient of a particular integral containing four- 
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point vertices. For some some solutions, the kinematics does not allow the coefficients of 
integrals with quartic vertices to be determined. For example, if the kinematic constraints 
due to three-point vertices force the spinors of two nearest neighboring legs of the four-point 
subamplitude to be proportional, Aj oc Xj, the sum of momenta of these legs will be on-shell, 
{ki + kjf = 0. This can place an internal propagator of the four-point subamplitude on-shell, 
where it diverges. This effectively selects out only those terms with this propagator present 
and loses the four-point contact contribution. In such cases, we determine the integral 
coefficient by using a different solution to the cut conditions. 

B. Evaluating the five- loop integral coefficients 

A maximal cut for determining the coefficient of integrals with a given set of propagators 
is of the form, 

n-4siT),. (7.1) 

h ^k=i ^ 

where h signifies the different helicity configurations and particle types that can contribute 
and m is the number of tree amplitudes appearing in the cut. In this equation the cut 
momenta li,l2, ■ ■ ■ ,lc are all on shell. Euler's formula relates the number of tree ampli- 
tudes that appear to the number of cut lines; at five loops m = c — 4. As discussed in 
section IVI Al a given kinematic solution to the cut conditions will allow only a subset of 
helicity configurations to contribute. 

The simplest situation is when an integral has only cubic vertices, as is true for the 
integrals of fig. [71 In this case one can always choose singlet kinematic solutions to the on- 
shell conditions so that only gluons propagate in each loop. After cutting all the propagators, 
one obtains the numerator subject to the cut conditions, 

^ = (n^8r) ' (/?,/2^•••,4=o) (7.2) 

Vj—]^ / singlet 

where the tree amphtudes ^(J)'^ are purely gluonic. In some cases, such as cuts isolating 
integrals Ji or I2, only a single term appears, but in others a sum of terms appear. For 
example, there are five contributions to the cut with propagator configuration of I21 and 
J22; the integral I21 appears four times in the cut as well as in the amplitude, but with the 
numerator factor permuted, while /22 appears one time. 
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(a) 



FIG. 17: (a) The simplest cut determining the coefficient of I37. (b) A particularly good choice 
of cut for determining the coefficient of I^g. All lines are cut and carry on-shell momenta. In all 
loops except those indicated by a dashed circle, only gluons propagate; the dashed circle indicates 
that all particle types can circulate. Other allowed helicity configurations are obtained by flip- 
ping all helicities in these loops. Grey blobs represent four-point amplitudes, which hide possible 
propagators inside. 

C. Examples of evaluations of integral coefficients 

As a first example of the determination of the coefficient of an integral, consider the 
maximal cut of the "ladder" integral Ji. Here = ist^A^^'^^{l, 2, 3, 4), which is independent 
of the loop momenta. This coefficient was determined long ago, from iterated two-particle 
cuts [23]. To confirm this result using maximal cuts we solve the on-shell constraints If of 
all sixteen cut propagators. For the external helicities (1^,2^,3^,4+) there are 166 distinct 
kinematic solutions of which 62 correspond to singlet configurations. For (1^,2^,3^,4"'") 
there are 270 solutions of which 27 are singlets, and for (1~, 2"^^, 3"^", 4~) the numbers are 
similar: 270 solutions, with 28 singlets. We have confirmed that all these singlet solutions 
individually agree with the known result, providing a non-trivial check. This determines that 
the pseudo-conformal integral Ji enters with an overall factor of —1/32 in the normalized 
amplitude M^^^ given in eq. (15. ip . after accounting for the normalization conventions of the 
integrals in eq. (15.21) . (The 1/32 prefactor in eq. (15.11) does not appear in the product of tree 
amplitudes making up the cut, but appears in the L-loop amplitude, due to our convention 



38 



of including a factor of 2^ in eq. fl2.1l) .) In the integrals of fig. [7] and fig. [S] we have not 
included the overall factor of —1/32, but leave it as an explicit overall factor in eq. (15. ip . In 
the remaining part of this section, we will refer only to signs relative to Ji. 

As a second example, consider /17 and lis in fig- as well as I^q in fig. [TOl They have 
the same propagators and hence can contribute to same maximal cut . For external helicities 
2~, 3"^, 4+) there are 335 kinematic solutions of which 62 are singlets. The helicity 
configuration 2+, 3~, 4+) has 339 solutions and 48 singlets whereas (1^, 2+, 3+, 4~) has 
304 solutions and 102 singlets. Again all singlet solutions individually give results consistent 
with the coefficient of I4Q vanishing and both lu and Jig entering the amplitude with a 
numerical coefficient of +1 relative to Ji. 

As a third example, consider a maximal cut of integral J21. Together with integrals having 
the same propagators, I22, I35 and J49, there are nine potential terms to the numerator 
after symmetrization. For this cut the helicity configuration (1^, 2^, 3+, 4+) has 376 solutions 
of which 98 are singlets and (1~, 2+, 3~, 4+) has 384 solutions of which 58 are singlets. Note 
that hehcity (1~, 2+, 3"*", 4~) is related to (1^, 2", 3"*", 4+), by symmetry. It turns out that 
the singlets can never be made consistent with numerators of type I^g, hence its coefficient 
must vanish. But unfortunately, the singlet cuts do not uniquely fix the coefficients of the 
remaining integrals. For example, if we were to assume relative numerical coefficients of 
±1, there are exactly two possibilities, one involves five terms given by symmetrizations of 
numerators of type I21 and I22 and the other involves two terms of type J35. To resolve this 
situation we must instead consider cuts with fewer cut conditions imposed, to reduce the 
degeneracy of the kinematics. 

Maximal cuts of diagrams involving non-cubic vertices are only a bit more complicated. 
Luckily almost all cuts needed to determine the coefficients of the integrals have singlets in 
their solution set. Only J37, 139, J55 and I^g in figs. [9] and [10] do not have singlet solutions. 
For these cases we must use non-singlet cuts, such as those in fig. [T71 

As an example of a singlet solution with a quartic vertex, consider a maxi- 
mal cut of J32. An expression that correctly matches the singlet maximal cut is 
C^^ = i A^^'^'^ {1,2, 3, A) {s'^t'^ + ...), where ". . ." stands for 14 rational terms obtained from 
integrals Iq, In, lu, hi, I22 and /31, which also contribute to this cut. Since the coefficients 
of these integrals can be determined from other cuts, we simply subtract their contributions 
allowing us to determine the coefficient of 132 to be +1. There are now fewer cubic vertices in 
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the cut and consequently the number of kinematic solutions also drops: The three inequiv- 
alent external helicity arrangements each have 18 solutions that are not degenerate in the 
two four-point blobs. However, they differ in their singlet content: Helicities 2~, S"*", 4"*") 
have four singlets, 2+, 3~, 4"^) have no singlets, and 2"*", 3^, 4~) have exactly one sin- 
glet. (When solving for the kinematics, we do not include solutions which do not allow us to 
determine I^2i because the four-point contact terms are lost, as mentioned in section IVII A[ ) 

As mentioned above, integrals 137, Jag, Jss and J59 have no singlet solutions in their 
maximal cuts. For J37 an useful choice is to force scalars and fermions to circulate in only 
two independent non-overlapping loops; there is only a single kinematic solution with this 
property, which helicity configuration shown in fig. [T7l(a). J39, J55 and J59 all have simpler 
solutions with only one loop that carries fermions and scalars. For J39, for example, three 
kinematic solutions exist with this property, one of which is displayed in fig. [T7l(b). 

A cut where a fermion or scalar can circulate in only one of the loops takes the form, 

cl-'-^=^' E {fn^o-r) -4(ff^gr) +3(ff^sr) | (7.3) 

I i=l /gluon i=l /fcrmion i=l /scalar-* 

where h is the helicity of the particle in the unique loop with fermions and scalars. (Helicity 
is conserved along this loop, which given our all-outgoing convention means that it flips 
going from one vertex to the next.) For the complex scalars, the two helicities correspond 
to particle and antiparticle. For the cut in fig. [TTT b) only four vertices involve particles 
other than gluons, hence five y4*'''^*^'s can be pulled out of the sum as a common factor. The 
required four-point tree amplitudes for different particles follow from the supersymmetry 
Ward identities 40|], which are described, for example, in Appendix E of ref. 24|. From this 
cut we find that 139 does not contribute to the amplitude. Likewise the maximal cuts of Jss 
and /59 shows that coefficients of these integrals vanish. 

For the cut in fig. [TW a) one can arrange the kinematics so that the two loops that can 
carry fermions or scalars do not intersect, simplifying their evaluation. The structure is 
similar to eq. (17.31) . except that there are two independent sums over fermions, scalars and 
gluons. We will not present it explicitly, but instead just give the kinematic solution needed 
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for this cut, 



h — P — XpXp , h — Q — XqXq , 

, _ , r {p + q + h + k2f _ ,r{P + q + h + ki 

'3 — —^n^n- ; ; ; — ; r , U — — A„A 



h — ~XvXq 



^ {q \p + ki + k2\q ) ' ^ (g | p + fci |r ) ' ^ 



\2 



{v^\p — ki \q~) 

Here p and q are arbitrary null vectors in four dimensions and and are spinors corre- 
sponding to arbitrary null vectors r, v. The remaining seven loop momenta can be obtained 
by momentum conservation. 

This cut is the least discriminating one needed for fixing the coefficients of the integrals 
in the five-loop amplitude, and hence it contains the most terms. There are 79 terms of 
the right conformal weight that are candidates for the left-hand-side of (17. ip . but of these 
only 28 terms contribute to the amplitude. These terms are obtained from integrals J5, Jig, 
I2O) ^21, I22 and /27; the coefficient of L^j must thus vanish. Interestingly, this is the only 
maximal cut where integrals (/21 and I22) enter twice compared to their appearance in the 
amplitude. 

In some cases the maximal cuts cannot distinguish between different integrals, due to the 
degenerate nature of the kinematics. As an example, consider the maximal cut of /21, /22 
and J35 described above, which has two possible numerator combinations satisfying the cut 
conditions. On the maximal cut of these diagrams we find, 

(/21 + I22 - I35) = . (7.5) 

cut 

The combination of I21 and /22 make one possible numerator choice and 135 is another 
consistent choice with this maximal cut. We have checked that more than 700 kinematic 
solutions of this cut fail to distinguish between the possibilities. To resolve this situation we 
use less degenerate kinematics with fewer cut conditions imposed. The two cuts in fig. [TTl 
for example, resolve this ambiguity. 

This type of ambiguity can even affect combinations of integrals with different sets of 
propagators. As a nontrivial example, the following combination of integrals vanishes in all 
maximal cuts we have evaluated, other than the ones in fig. [17] and the maximal cut of /50: 

(/21 + /22 - hr + /31 + h2 + h3 + hi - hb - he + Iss) ^^^= . (7.6) 

This equation as well as eq. (17.51) should be interpreted as a recipe for determining a com- 
bination of terms that can vanish in a maximal cut; if a cut picks up any integral or its 
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permutation it should be included. (Note that the signs shown in figs. [7] and [S] are included 
in the definition of these integrals.) Note that eq. (17. 5p is the same ambiguity as eq. (17.61) . 
but restricted to cuts of /21's topology. Other than (17. 6p we have found no ambiguity that 
holds for all singlet solutions of a maximal cut. In any case, it can be resolved by using 
cuts with fewer on-shell conditions. In particular, the cuts in fig. [17] resolve the ambiguity 
(17. 6p . These cuts are only consistent with integrals /21, 122^ hj, hi: h2, and I34 included 
in the amplitude and I^q and J38 excluded. It is likely that this kind of ambiguity also 
exists at all higher loops when using maximal cuts, but again it should be resolved by using 
less-restrictive cuts. 

Although it is straightforward to solve analytically for any given kinematic configuration, 
it can get quite tedious since many of the pseudo-conformal integrals have well over 100 
singlet cuts each. It is therefore simpler to do so numerically. The bispinor formalism, 
which is automatically on-shell, enables us to choose which kinematic solution to solve for 
numerically : Our procedure is to first assign spinors, one to each three-point vertex, with 
A's assigned to each vertex and A's assigned to each blob. If two nearest-neighbor 
three-point vertices are of the type, the two A spinors of the vertices are set equal to 
each other. The on-shell constraints force the A spinors to be proportional to each other, 
but we use the rescaling freedom of the spinors, mentioned below eq. (16.101) . to set the 
proportionality constant to unity. Similarly, if two nearest-neighbor three-point vertices are 
of the type, the A spinors are set equal to each other. This gives us a list of A's and A's 
that uniquely specifies the solution, but their values are not yet determined. The momentum 
of a cut propagator between blobs of opposite sign is given by taking the tensor product of 
the two spinors associated with the blobs the propagator connects to. One must also allow 
for a complex scale factor multiplying each momentum between these blobs since it is not 
possible to simultaneously remove the proportionality constants in both A and A. We solve 
the momentum conservation constraints by numerically minimizing the sum of the squares 
of absolute value of all momentum conservation relations that should vanish. At the solution 
point this vanishes. 

In some cases the numerical convergence is insufficiently fast. If necessary a given cut can 
always be analyzed analytically. But it is simplest to discard unstable or poorly convergent 
solutions because there are plenty of other solutions available. We have performed on the 
order of 100 singlet maximal cuts corresponding to each of the propagator configurations of 
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FIG. 18: The two singlet cuts containing only gluons in D = 4. 



the pseudo-conformal basis integrals with only three-point vertices, effectively exhausting 
the singlets. For the integrals also containing four-point vertices, we have performed on the 
order of 10 singlet maximal cuts, again effectively exhausting the singlets. We have also 
checked a handful of non-singlet solutions, in particular for those diagrams which have no 
singlet solutions. In all cases, we find that the cuts are consistent with the ansatz (15.11) . 

VIII. CROSS CHECKS ON COEFFICIENTS FROM TWO- AND THREE- 
PARTICLE GENERALIZED CUTS 

The kinematics used in the maximal cuts is rather restricted, so additional checks are 
desirable. We have evaluated two such generalized cuts in four dimensions. In D dimensions 
we evaluated various two-particle cuts. These cuts also provide a confirmation that no other 
integrals appear in the amplitudes besides the pseudo-conformal ones. 

A. Cuts in four dimension 

The easiest four-dimensional cuts to evaluate are the singlet cuts, involving only MHV 
gluon tree amplitudes and a single three-particle cut, shown in fig. [131 As described in 
section IVll the non-75 terms, obtained after dividing by the tree amplitude (see eq. (16.161) ) 
give precisely half the cut of the final amplitude at least through four loops. By evaluating 
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the cut of fig. [TSl we obtain a non-trivial check, since we find that a similar result holds at 
five loops. 

Our evaluation confirms agreement of the non-75 terms in the singlet cut of fig. [TST a) 
with 1/2 the value of the unintegrated corresponding cut of the ansatz (15. ip . This provides 
a non-trivial confirmation that we have properly determined the coefficient of integrals 

-^3; -^8; -^9) ho, -^13) hi, hli hs, and/28, (8.1) 

from the maximal cuts. Similarly, we have confirmed that the non-75 terms in the singlet 
cut (b) agrees with 1/2 times the value of the corresponding cut of the ansatz (15.11) . This 
checks that the coefficients of the integrals, 

hi hi hi ■ ■ ■ 1 h2i hii hsi hbi 

and J30 , (8.2) 

are also correct. Moreover this also checks that integrals which have cuts of the forms in 
fig- [HI but are not pseudo-conformal do not appear in the amplitude. 



B. Cuts in D dimensions 



A more rigorous check comes the evaluation of the D-dimensional cuts. As already 
mentioned, beyond one loop, no theorem has been proven that four-dimension cuts are 
sufficient for determining complete amplitudes in supersymmetric theories. It is therefore 
important to evaluate at least some cuts in D dimensions. This is especially true if we wish 
to apply the results away from D = A. 

We evaluate the D-dimensional cuts of MSYM, by interpreting it instead as ten dimen- 
sional M = 1 supersymmetric Yang-Mills, dimensionally reduced to D dimensions. As 
mentioned in section IVIl this way of evaluating the MSYM amplitudes has the advantage 
of simplifying the bookkeeping on which states are present: the M = 1 multiplet consists of 
only a single gluon and gluino, each of which is composed of 8A''c degrees of freedom. With 
this formulation, all states are included, using Z)-dimensional momenta in the cuts. 

The simplest class of integrals to check in D dimensions are ones which can be constructed 



by iterating two-particle cuts, following the discussion of refs. 23|, \2S\- The two-particle 
sewing equation, which is valid in D dimensions, is. 



Ar(^i,i,2,/2)Ar(-/2,3,4,-/o = -^A*r(i,2,3,4)- — — -_. (8.3) 

N=A states I'l ~ '^^) I'l + '^4j 
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(a) (b) (c) 

FIG. 19: The D-dimensional two-particle cut dividing the five-loop amphtude into (a) two two-loop 
amplitudes, (b) a one-loop and three-loop amphtude and (c) a four-loop and a tree amphtude. Ah 
physical states are summed over in the cuts. 

Since the tree amplitude A^^^'^ appears on the right-hand-side, the same two-particle sewing 
algebra appears at the next loop order. The iterated two-particle cuts allow us to confirm 
that the coefficients of integrals 

Ji . . . Jio, /i5, /i6, ^19, ho, hi ■ ■ ■ hn, hi, hs, hi, hs , (8.4) 

have all been determined correctly. 

We have also checked the D-dimensional two-particle cuts which split the five-loop four- 
point amplitude into a product of a two two-loop amplitudes, a one-loop and three-loop 
amplitude and a four-loop and a tree amplitude, as depicted in fig. [191 Using D-dimensional 
cuts we have evaluated the coefficients of all integrals appearing in two- and three-loop 
amplitudes, leaving the external legs in D dimensions. These are all proportional to the 
D-dimensional tree amplitude. We may likewise use the D-dimensional four-loop amplitude 



subject to the same assumptions make in ref. [10(], namely the absence of certain triangle 
subintegrals and the appearance of the tree-level kinematic tensor as an overall coefficient. 
We can then apply the two-particle cut sewing equation (18. 3p to con&m the coefficients of 
various five-loop integrals. This allows us to provide additional checks via D-dimensional 
unitarity that integrals, 

h ■ ■ ■ ho, hb ■ ■ ■ ho, hs, hb, he, hi ■ ■ ■ hr, hi ■ ■ ■ hs , (8-5) 

all have the coefficients presented in section I VI II 

To have a complete proof that the ansatz (15. ip is complete, one would need to to conffim 
from D-dimensional unitarity that these remaining integrals enter with the coefficients de- 
termined in section IVIII and that there are no other (non-conformal) integrals present. We 
leave this for future work. 
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In general, it is rather surprising that four-dimensional unitarity cuts are sufficient to 
determine the amplitudes in all dimensions. The maximally supersymmetric theory, how- 
ever, is special. Our D-dimensional study here provides non-trivial evidence that at least 
at four points, the four- dimensional cuts suffice. This result may be understood as a direct 
consequence of only pseudo-conformal integrals being present, with coefficients independent 
of the number of dimensions. 



IX. CONCLUSIONS 



In this paper we presented an ansatz for the five-loop four-point planar amplitude of 
maximall y s upersymmetric Yang-Mills amplitudes in terms of a set of pseudo-conformal 



integral 
gators 



26 



10 



291] . We introduced a method based on cutting the maximal number of propa- 



in each integral, to determine very efficiently the coefficients of the integrals 



22 1 with less restrictive 



as they appear in the amplitude. We then used generalized unitarity 
cuts, both in four and D dimensions, to verify the correctness of the expressions determined 
in this way. 

Our ansatz for the planar five-loop four-point amplitude relies on a basis of pseudo- 
conformal integrals, and assumes that the amplitude can be expressed entirely in terms 



of such integrals. T 
conformal integrals 



lese integrals are the dimensionally-regulated counterparts of off-shell 



10 



291, limited to those which have logarithmically-divergent on-shell 



limits. This assumption has been tested and conffimed by explicit calculation through four 
loops. The assumption provides a compact basis of (plausibly independent) integrals, and 
reduces the problem of computing the amplitude to that of determining the coefficients of 
each integral. We have provided strong evidence that this continues to hold through at least 
five loops, through the evaluation of a large variety of generalized unitarity cuts, including 
ones evaluated in D dimensions. The computation of additional cuts in D dimensions would 
make it possible to prove that our expression is indeed complete. Alternatively, we may 
wonder whether it is possible to link the conformal invariance of the theory to the absence 



of non-conformal integrals. An important cross-check would come 
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Tom showing that the 



infrared singularities of the amplitude have the predicted form 

The set of integrals that appears in the expression for an L-loop MSYM amplitude is a 
subset of all pseudo-conformal integrals. It is interesting that the integrals which do appear, 
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do so with coefficients ±1. We presented heuristic rules which give a partial understanding 
of the signs of these coefficients. It would of course be very useful to have a complete set of 
heuristic rules for predicting all signs and zeroes to arbitrary loop order. 

This maximal form of generalized unitarity we employed should also prove useful for 
determining non-planar contributions. For example, the non-planar contributions to the 
subleading-color three-loop amplitude shown in fig. [2D] are easily determined from maximal 



cuts. These contributions are in agreement with known results [2J, |3J|. (Note that with the 
cut conditions imposed (/ + k^)'^ and 21 ■ are indistinguishable in the second integral of 
fig- [20l so other cuts are necessary to determine the proper factor.) 

Our determination of coefficients also relied on special properties of the four-point ampli- 
tude. How can one compute amplitudes with a larger number of external legs? While some 
extension to the techniques presented in the present paper will certainly be necessary, they 
provide a very good starting point. In the planar two-loop five-point amplitude j^, U^, for 
example, terms with even parity relative to the tree amplitude also appear to be expressible 
purely in terms of pseudo-conformal integrals. The parity-odd terms require further study. 

Beyond computations of gauge-theory amplitudes, the maximal-cut method described 
here should also be useful in higher-loop studies of quantum gravity. Recent calculations have 
established 3J] that the three-loop degree of divergence in four dimensions (or equivalently 
the critical dimension) of A/" = 8 supergravity is — contrary to widely-held expectations 
— the same as that of A/" = 4 supersymmetric gauge theory. There are other indications 



that the supergravity theory may even be ultraviolet finite beyond three loops 
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33|, 



36|. These investigations point to the need for higher- loop computations of supergravity 



amplitudes, in order to establish the critical dimension in which they first become ultraviolet 
divergent. In the approach advocated in ref. (24], cuts of MSYM gauge-theory amplitudes 
can be used to construct cuts of A/" = 8 supergravity amplitudes. The present paper provides 
the required planar amplitudes at five loops. The non-planar contributions are more difficult, 
but should be within reach. This task would be considerably easier if a non-planar analog 
of the pseudo-conformal integrals were identified. 

Our expression for the planar five-loop four-point MSYM amplitude presented in this 
paper has two obvious applications. The first would be the extraction of the planar five- 
loop cusp anomalous dimension, allowing a further check of the conjectures of refs. 0, lo| . 
Another application would be a five-loop check of the iterative structure of the amplitude. 
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FIG. 20: Non-planar examples of three- loop integrals confirmed by cutting all the propagators. 
These agree with the results of ref. 24 ] . 



This would provide a rather strong check of the all-loop resummation of maximally helicity 
violating amplitudes proposed in refs. js, 4|, and help reinforce a link to a recent string-side 
computation of gluon amplitudes 12(]. The latter computation, together with all- loop-order 
resummations, opens a fresh venue for quantitative studies of the AdS / CFT correspondence. 
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